Symmetry-dependent phonon renormalization in monolayer M0S2 transistor 
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Strong electron-phonon interaction which limits electronic mobility of semiconductors can also 
have significant effects on phonon frequencies. The latter is the key to the use of Raman spec- 
troscopy for nondestructive characterization of doping in graphene-based devices. Using in-situ 
Raman scattering from single layer M0S2 electrochemically top-gated field effect transistor (FET), 
we show softening and broadening of Ai g phonon with electron doping whereas the other Raman 
active mode remains essentially inert. Confirming these results with first-principles density 
functional theory based calculations, we use group theoretical arguments to explain why Ai g mode 
specifically exhibits a strong sensitivity to electron doping. Our work opens up the use of Raman 
spectroscopy in probing the level of doping in single layer MoS2-based FETs, which have a high 
on-off ratio and are of enormous technological significance. 

PACS numbers: 78.30.-j 



Discovery of graphenei stimulated an intense research 
activity due to interesting fundamental phenomena it ex- 
hibits as well as the techonological promise it holds in a 
broad range of applications ranging from sensors to nano- 
electronics. Vanishing bandgap of a single layer graphene 
is a sort of a limitation in developing a graphene-based 
field effect transistor with a high on/off ratio. This has 
spurred efforts to modify graphene to open up a gap and 
towards development of other two dimensional materials 
like M0S2, WS2 and boron nitride (BN), both experimen- 
tally and theoretically. Avenues to open up gap through 
modification of graphene include quantum confinement in 
nanoribbons 2 , surface functionalization 3 , applying elec- 
tric field in the bilayer— deposition of graphene on other 
substrates like BN 6,7 , and B or N substitutional doping 8 , 
which require fine control over the procedure of synthesis. 

In contast to graphene, single layer M0S2 consisting of 
a hexagonal planar lattice of Mo atoms sandwiched be- 
tween two similar lattices of S atoms (S-Mo-S structure) 
with intralayer covalent bonding is a semiconductor with 
a direct band gap of ~ 1.8 eV, and is quite promising 
for FET devices with a high on-off ratio. It has been 
shown that the luminescence quantum yield of monolayer 
M0S2 is higher than its bulk counterpart^^. Recently a 
monolayer M0S2 transistor 11 has been shown to exhibit 
an on-off ratio of ~10 8 and electron mobility of ~200 
cm 2 /V-sec. These values are comparable to silicon based 
devices and make M0S2 based devices worth exploring 
further. It is known that in a field effect transistor, car- 
rier mobility is limited by scattering from phonons and 
the maximum current is controlled by hot phonons. Both 
these issues in a FET depend on the electron-phonon cou- 
pling (EPC). Raman spectroscopy has been very effective 
to probe EPC for singlei^r— and bilayer graphene — ~— 
transistors by investigating the renormalization of the G 
and 2D modes as a function of carrier density. 

Recent layer-dependent Raman studies of single and 
few layers of M0S2 18 have shown that freuqency of E^ 
phonon increases as the number of layers decreases, 
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FIG. 1: (Color online) (a) Schematic of experimental set up. 
Source, drain and gate electrodes were marked as S, D and G 
respectively. The dotted box is the polymer electrolyte layer, 
(b) AFM height profile of the monolayer M0S2 flake. (Right 
Inset) AFM image showing the dotted line along which the 
height profile is taken. Scale bar is 1 fim. (Left Inset) Atomic 
displacements, indicated by arrows, corresponding to E^ and 
Ai g modes as viewed along the [1000] direction, (c) Ids-Vds 
plot at various top gate voltage Vtg- (d) Ids as a function 
of Vtg at Yds = 200 mV. The on-off ratio was - 10 5 and 
the device mobility was estimated to be 50cm 2 /V-sec. Inset 
shows the transfer characteristic in linear scale. 



whereas the frequency of Ai g phonon decreases 18 . This 
has been recently explained in terms of enhanced di- 
electric screening of the long range Columb interaction 
between the effective charges with increasing number of 
layers^. For non-destructive characterization of carriers 
in recently developed mono-layer M0S2 transistor—, Ra- 
man spectroscopy can be quite useful and requires precise 
knowledge and understanding of phonon renormalization 
of single M0S2 layer as a function of carrier concentra- 
tion. The temperature dependence of the mobility in n- 
type bulk M0S2 had been attributed to the scattering of 
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carriers by optical phonons that modulate thickness 20 im- 
plying that the Ai g phonons with atomic displacements 
parallel to the c-axis should be involved in controlling 
mobility of the carriers. 

In this letter, we report in-situ carrier dependent Ra- 
man study of a top gated single layer M0S2 transistor 
achieving a maximum electron doping of ~2x 10 13 /cm 2 . 
The transfer characteristic of the top gated device shows 
an on-off ratio ~ 10 5 and a field effect mobility of 50 
cm 2 /V-sec. We show that the Ai g mode shows a strong 
doping dependence; the phonon frequency decreases by 4 
cm -1 and linewidth broadens by 6 cm -1 for electron dop- 
ing of 1.8 xl0 13 /cm 2 . Phonon frequency and linewidth 
of the El g mode show a much less dependence on car- 
rier concentration. The difference in the behaviour of 
these optical phonons are explained quantitatively using 
density functional theory (DFT). 

Fig. QJi shows a schematic of our experimental set 
up. Single layer M0S2 flakes were mechanically exfoli- 
ated from a bulk single crystal procured from M/s. SPI 
Supplies and transferred on a 300 nm Si02 on a degener- 
ately doped p-type silicon substrate (procured from M/s. 
XT Wafer). After optical identification, the flake height 
is measured by contact mode AFM to be 0.7 nm, in 
agreement with the S-Mo-S layer thickness (see Fig.[T|3). 
Standard electron-beam lithography and deposition were 
done to form ~50 nm thick Au contacts as source (S) , 
drain (D) and gate (G) electrodes. Room temperature 
Raman spectra were recorded with 514.5 nm laser exci- 
tation with Witec confocal spectrometer using 50X long 
working distance objective. Laser power was kept below 
lmW to avoid sample heating. Electrical measurements 
were done with Keithley 2400 source meters. For top 
gating, we have used solid polymer electrolyte compris- 
ing of a mixture of LiClCU and polyethylene oxide (PEO) 
in the weight ratio 1:8. Solid polymer electrolyte as the 
gate material was chosen because of its high capacitance- 
enabling high carrier concentration with low gate voltage 
(< 2V). At the same time, being almost transparent, it 
allows us to perform in-situ optical measurements with 
simultaneous electrical characterization. Fig. shows 
drain-source current (Ids) as a function of drain-source 
bias (Vds) f° r a representative device of length (L) 2.5 
/xm and width (W) 1.5 /am. A small non linearity may be 
due to the Schottky barrier at the contacts. The transfer 
characteristics of the device in Fig.[T]i is plotted in a semi 
log scale. The on-off ratio for our device is ~10 5 and the 
low field effect mobility (fi) of our device was calculated 
(/i = 1^7 v d 9 sCtg ) to ^ e ^ cm 2 /V-sec. Here g m is the 
transconductance (g m = |^P^ ) of the device. Enhanced 
photoconductivity was observed when the incident laser 
beam was focused on the sample at gate voltages of 0.5, 
1.0, 1.5 and 2.0V in agreement with recent reports^. 
The photoconductive response was large at lower gate 
voltages but decreases at higher gate voltages when the 
conduction band starts getting populated. Raman spec- 
trum at each gate voltage was recorded only after the 
stabilization of the channel current. The threshold gate 



voltage Vt, at which the device switches from 'off' to 
'on' state was ~0.1 V. The gate induced electron concen- 
tration n is estimated using ne = Ctg(Vtg — Vt)- The 
value of Ctg is 1.5 /iF/cm 2 . In brief, we extracted the 
value of top gate capacitance Ctg from an experiment 
with bilayer graphene transistor involving dual gate con- 
figuration—. The back gate dielectric was 300 nm Si02- 
Several back gate sweeps were done for fixed top gate 
voltages. Since a maximum in the resistance corresponds 
to the charge neutrality point, the carrier concentration 
induced by both the gate voltages were equated to get 
the Ctg to be 1.5 /iF/cm 2 . Since the electrolyte in the 
present experiments is same as in Ref.[5], we take this 
value in estimating n. 

Fig. [2^ shows the evolution of zone center phone E^ 
and A\ g modes of M0S2 monolayer at different top gate 
voltages. As depicted in Fig. [Tb inset, A\ g phonon in- 
volves the sulphur atomic vibration in opposite direc- 
tion along the c axis (perpendicular to the basal plane) 
whereas for E\ g mode, displacement of Mo and sulphur 
atoms are in the basal plane. Lineshape parameters 
were obtained by fitting a sum of two Lorentzians to the 
data. Fig. [2J) and Fig. [2fc show the shift of the mode fre- 
quencies and the corresponding full width at half maxi- 
mum (FWHM), respectively, as a function of gate volt- 
age. The dependence of the change in mode frequencies 
[Auj = u(n ^ 0) — uo(n = 0)] on the carrier concentration 
(n) is shown in Fig. [2]i. For a maximum electron concen- 
tration of 1.8 x 10 13 /cm 2 , the Ai g mode frequency soft- 
ens by 4 cm -1 , as compared to only ~ 0.6 cm -1 for the 
El g mode. The linewidth of the Ai g mode increases sig- 
nificantly by ~6 cm -1 for the maximum doping achieved, 
whereas the linewidth of the E^ mode does not show 
any appreciable change. These results show that the A± g 
phonon renormalization could be used as in-situ read-out 
of the carrier concentration in M0S2 devices. We will now 
quantitatively understand the different renormalization 
of the two modes, A± g and E^ due to electron-phonon 
interaction. 

Our calculations are based on first-principles density 
functional theory as implemented Quantum Espresso 
package, 21 a plane- wave basis set (70 Ry cutoff) 
and norm conserving pseudopotentials (Rappe-Rabe- 
Kaxiras-Joannopoulos 22 ). Exchange-correlation energy 
of electrons is approximated with a local density approx- 
imation and a parametrized functional of Perdew and 
Zunger 23 . A monolayered form of M0S2 is simulated us- 
ing a periodic supercell, with a vaccum of - 15 A sep- 
arating adjacent periodic images along the z direction. 
Integrations over the Brillouin zone were sampled with 
uniform 24 x 24 x 1 and 48 x 48 x 1 k— point meshes in 
calculation of total energy and electron-phonon coupling, 
respectively. 

Our calculations show that monolayer M0S2 is a di- 
rect band gap semiconductor with a band gap of 1.8 eV 
at K-point, in good agreement with known results 24 . Va- 
lence and conduction band edges primarily consists of Mo 
Ad states with some hybridization with S 3p states [see 
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FIG. 2: (Color online) (a) Raman spectra of monolayer M0S2 at different top gate voltages Vtg- Open circles are experimental 
data points, the gray (red) lines are Lorentzian fits to the total spectrum and the black color lines are the Lorentzian fit to 
individual peak. Change in the (b) phonon frequency Acj and (c) FWHM of A\ g and E^ g modes as a function of Vtg- Change 
in zone center phonons Auj of (d) from experiment and (e) from DFT calculations as a function of electron concentration n. 
(f) Electron-phonon coupling of Ai^ and modes as a function of n. (g) Ratio of EPC [^A lg {ft 7^ 0)/Aa 19 {tl — 0)] shown by 
open squares and phonon linewidth [FWHMa 19 (n 7^ 0)/FWHMa 19 (n = 0)] shown by filled squares as a function of n. 
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very good agreement with our experimental results shown 
in Fig. [2]i. In order to understand this trend, we carry 
out a systematic study of electron-phonon coupling as 
a function of electron doping. The electron-phonon cou- 
pling (EPC) of a mode v at momentum q (with frequency 



is calculated as 25 



FIG. 3: (Color online) Contour plots (in light gray) of charge 
density (\ip(r)\ 2 ) (a) at the top of the valence, and (b) bottom 
of the conduction bands at K point. While states near the 
top of valence band consist of Mo 4d states with some hy- 
bridization with S 3p states, those near the conduction band 
edge have clearly the character of d z 2 state of Mo. The dark 
gray (red) and black (blue) spheres are Mo and sulphur atoms 
respectively. 
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where uj and N(ef) is the phonon frequency and elec- 
tronic density of states at the Fermi energy, respectively. 
The electron-phonon coupling matrix element is given by 



Fig. [3] a and b]. In particular, states near the bottom 
of conduction band near the K-point have a character 
of d z 2 state of Mo, and that at the top of valenec band 
has character. We note that charge density (\ip(r)\ 2 ) 
associated with each of these states has a full symme- 
try of M0S2 layer, which has important consequences for 
electron-phonon coupling (to be elaborated later). 

We simulated electron doping in M0S2 by adding a 
small fraction of electrons to its unit cell. Doping has 
contrasting effects on the frequencies of Ai g and E^ 
optic modes as shown in Fig. [2^. While the former 
mode is found to soften significantly (~ 7 cm -1 at 
~ 1.8 x 10 13 /cm 2 doping), the latter is hardly affected, in 



5Z + «' fc = (2^) 2 <^ + ^l A VI«, (2) 

where iph.j is the electronic wavefunction with wavevec- 
tor k and energy eigenvalue e^j for band j and M is the 
ionic mass. AV qu is the change in the self-consistent po- 
tential associated with a phonon of wavevector qr, branch 
v and frequency uo qv . Eq [2] defines the scattering of an 
electron from band j to band i due to the phonon v with 
momentum q. Our results in Fig. [2f show that Ai g mode 
couples much more strongly with electrons than the E^ 
mode. This can be understood using group theoretical 
analysis of symmetry. Ai g mode has a symmetry of the 
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lattice (the identity representation, i.e. the structural 
distortions in this mode do not break the symmetry of 
M0S2, see Fig.[T]3 inset). As a result, all electronic states 
can have a nonzero expectation value in Eq[2]for the per- 
turbation of Ai g mode, giving a large electron-phonon 
coupling in Eq [TJ Electron doping leads to occupation 
of the bottom of the conduction band at K-point states 
which have a character of d z 2 of Mo [see Fig. [3b]. The 
|^(r)| 2 of the states near K— point also transform accord- 
ing to the identity representation A\ g . Hence, changes 
in occupation of these states with electron doping yield a 
significant change in the EPC of the K\ g phonon. In con- 
trast, the matrix element in Eq [2] vanishes for E^ mode 
(orthogonality of A\ g and E^ representations) and its 
coupling with electrons is weakly dependent on doping. 
Frey et ali2£ have put forward similar conclusions. 

It is interesting to compare EPC of M0S2 with that of 
graphene. First of all, electron doping of about 1.8xl0 13 
cm -2 results in hardening of G-band of graphene — by 
about 10 cm -1 and softening of K\ g mode of M0S2 by 
about 4 cm -1 . The G-phonon renormalization occurs 
due to phonon induced electron-hole (e-h) pair creations. 
For the G mode (q ~ 0) in graphene, it involves e-h cre- 
ation within a valley. Doping the graphene blocks the 
generation of phonon induced e-h pairs and hence af- 
fecting the phonon self energy. Secondly, distortions of 
the structure with atomic displacements of a G-phonon 
lead to a mere shift of the centre of the Dirac cone (and 
in the Fermi-surface). If electrons follow these distor- 
tions remaining in their ground state (adiabatic limit), 
their energy cost does not change with doping and hence 
doping would not result in any shift of frequency of the 
G-band. However, it is the breakdown of adiabatic ap- 
proximation (the fact that electrons do not follow nu- 
clear motion remaining in their ground state) that is re- 
sponsible for the energy cost and hardening of the G- 
band with electron and hole doping — . Another way 
to explain is in terms of Kohn anomaly in graphene for 
q = phonon which is weakened on doping. In con- 
trast, in semiconducting M0S2, phonon renormalization 
occurs within the adiabatic approximation. The elec- 
tron doping results in occupation of the anti-bonding 
states in the conduction band of M0S2 making the bonds 



weaker and the K\ g mode, which preserves the sym- 
metry of the lattice, softens. While the EPC of Ra- 
man mode at Y point of M0S2 exhibits a strong depen- 
dence on doping, similar dependence is seen for K— point 
phonons of graphene 25 and our symmetry-based expla- 
nation applied here as well. We note that the latter 
transforms according to identity representation of the 
symmetry group at (or near) the K— point. The EPC 
argument applies to the linewidth as well. Two differ- 
ent mechanism contribute to the phonon linewidth: (a) 
the EPC contribution (FWHM^ PC ) and (b) the contri- 
bution (FWHM an ) arising from anharmonic effectsSLSS. 
We can express, FWHM = FWHM^ PC + FWHM an . 
The phonon linewidth, FWHM^^, is proportional to 
EPC associated with a particular mode 27 i 28 . The in- 
crease in Aip linewidth is a result of the strengthening 
of electron-phonon coupling (A) with doping (Fig. [2k and 
Fig. EF). Instead of the absolute values, we compare the 
ratios of FWHM (FWHM n ^ /FWHM n=0 ) and EPC val- 
ues (A n ^o/A n =o) for the K\ g mode (see Fig. [2^). The 
FWHM ratio follows the same trend as the EPC ratio 
establishing that the increase in linewidth is due to in- 
crease in the electron-phonon coupling values (A) with 
doping. 

In summary, we have demonstrated that electron dop- 
ing in single layer M0S2 results in softening specifically 
of its Raman active K\ g phonon, accompanied by an in- 
crease in the line-width of its Raman peak. In compari- 
son, the other Raman mode with E^ symmetry is quite 
insensitive to electron doping. This is due to a stronger 
electron-phonon coupling of the K\ g mode than of the 
Er>£ mode, confirmed with first-principles DFT calcula- 
tions and symmetry arguments. Our work shows how 
Raman scattering can be effectively used to characterize 
level of doping in gated FET device based on single-layer 
M0S2 with a high on-off ratio of ~ 10 5 and having po- 
tential to be used in digital electronics and sensors. 
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